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We make a first demonstration of a sub-keV electron recoil energy threshold in a dual-phase liquid 
argon time-projection chamber. This is an important step in a program to build a detector capable of 
identifying the ionization signal resulting from nuclear recoils at a few keV and below. We obtained 
this result by observing the peaks in the energy spectrum at 2.82 keV and 0.27 keV, following the K- 
and L-shell electron capture decay of 3r Ar. We describe the details of the 37 At source preparation, 
as this calibration technique may prove useful, e.g. for dark matter direct detection experiments. A 
55 Fe internal x-ray source was also measured simultaneously and provided another calibration point 
at 5.9 keV. We discuss the ionization yield and electron recombination in liquid argon at the three 
calibration energies. 



I. INTRODUCTION 

Dual-phase noble liquids detectors have become a pop- 
ular choice for WIMP dark matter experiments thanks 
to their target mass scalability, low background, and 
low threshold. More recently, they have been proposed 
also to search for coherent neutrino-nucleus scattering 
(CNNS). In both cases, the only signature in the detec- 
tor is a low- energy nuclear recoil. For the most common 
WIMP scenarios, nuclear recoil energies between 10-100 
keVr in liquid Ar or Xe are expected [T]. For CNNS, de- 
tection sensitivity to nuclear recoils of a few hundred eVr 
is needed, especially if reactor neutrinos are employed 
[2J. Energy thresholds of 5-10 keVr are commonly ob- 
tained in dark matter experiments with liquid xenon or 
argon. The XENON10 experiment recently showed that 
an energy threshold as low as 1 keVr, corresponding to 
only a handful of ionized electrons, might be obtained [3] 
in liquid xenon. This result was obtained from extrap- 
olation of the model described in [4]. It is essential to 
measure the low-energy response of noble gas detectors 
to allow their use for CNNS, and to probe broader ranges 
of light-mass WIMPs. 

Here we report the first demonstration of sub-keV 
electron-recoil spectroscopy in a dual-phase argon detec- 
tor. This is achieved by detecting the proportional scin- 
tillation produced in the gas phase by ionization electrons 
after they have been extracted from the liquid. Specif- 
ically, we observed the 270 eV cascade following L-shell 
electron capture in 37 Ar. We use the single electron sig- 
nal as an absolute calibration for the ionization chan- 
nel. We then suggest that low energy electron recoils can 
be modeled within an existing framework. The results 
shown here are a first step toward providing detector ca- 
pability to measure the ionization yield of nuclear recoils 
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at few keVr and below, as well as for enhancing the reach 
of axion search via the axio-electric effect. 

In liquified noble gas detectors, the energy transferred 
to the medium by a particle results results in ionization, 
scintillation and heat [5] . The fraction lost to each chan- 
nel depends on both the incident particle type and en- 
ergy. Indeed electron recombinations varies with ioniza- 
tion density along the track as well as the electric field. A 
discussion of the mechanism and yields of energy transfer 
in liquefied noble gases lies outside the scope of this pa- 
per. Extensive literature is available on the subject (see 
for example [6] and references therein). However data 
are still scarce at low energies for both electron and nu- 
clear recoils. In addition, the measurements have mostly 
focused so far on scintillation yields alone, rather than 
ionization. 



II. DETECTOR AND CALIBRATION SOURCES 

We have built a prototype dual-phase liquid argon 
time-projection chamber in the context of our program 
to detect coherent neutrino-nucleus scattering j7J. Our 
detector emphasizes high-sensitivity measurement of the 
ionization signal by means of secondary scintillation in 
the gas region. A particle interaction in the liquid phase 
produces primary scintillation (SI) and ionization [5]. 
The electrons are drifted away from the interaction site 
by an electric field and extracted into the gas where they 
create secondary scintillation (S2) that we detect, light 
collection was optimized for the S2 rather than SI signal, 
and as a result limits the detection of SI below about 10 
keV. 

A schematic view of the detector is provided in Fig. 
[T] It consists of ^100 g of liquid argon as the active vol- 
ume surrounded on the bottom and sides by ~1 kg of 
inactive liquid. A gaseous region extends on top of the 
liquid and is kept at a constant pressure of ^820 torr. 
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FIG. 1. Schematic drawing of the vertical cross section of the 
detector central region. The four copper rings are negatively 
biased and shape the electric field in the drift region. The 
anode is at ground potential. A SS mesh (50 mesh/inch, 
30 /im dia wire, 88% optical transparency) is used for the 
anode and cathode grids. The extraction grid is obtained 
using 15 /xm SS wire with 1 mm spacing. The scintillation 
light is collected by four 1" PMTs through a TPB-coated 
silica window. The electric fields used to obtain the data 
reported here are also shown. 



The temperature of the gas is that of the Ar vapor at 
that pressure. The active volume is 3.65 cm high, and 
5 cm in diameter. A set of four copper rings delimits 
this active volume and shapes the electric field needed 
to drift the electrons in the liquid. An extraction grid is 
placed 5-10 mm below the liquid level, while the anode is 
located 2.5 cm above the extraction grid. Electric fields 
up to 4 kV/cm in the drift region and up to 10 kV/cm 
in the gain were achieved. The argon scintillation light 
is detected through a fused silica window by a squared 
array of four 1" Hamamatsu R8520 photomultiplier tube 
(PMT) modified for cryogenic operation. A 0.05 mg/cm 2 
coating of TetraPhenyl Butadiene (TPB) over the silica 
window acts as wavelength shifter for the Ar UV scintil- 
lation light [S] . The PMTs signals are triggered and digi- 
tized at 500 MHz using a LeCroy Waverunner 104MXi-A 
8 bit digital oscilloscope. 

The detector response was studied using different 
gamma sources. A 500/iCi 241 Am gamma source was col- 
limated to produce a thin (~2 mm) beam spanning the 
entire drift region of the detector. Photoelectric events 
(60 keV) from the 241 Am were used to assess argon pu- 
rity by measuring the electron lifetime extracted from 
the dependance of the total ionization signal from event 
depth in the liquid. The latter is calculated using the 
time delay between the SI and S2 scintillation signals. 
For the measurements reported here, no appreciable loss 
of electrons due to impurities was observed within the 
detector's energy resolution for events produced at dif- 



ferent depths. The 241 Am purity data were taken at a 
drift field of 3.0 kV/cm. 

Two sources located within the active liquid volume 
were used for calibration of the detector to low energy 
electron recoils. A 5 mm diameter electroplated 55 Fe bare 
source provided 5.90 keV and 6.49 keV x-rays with ~ 1 
kBq activity. The source was mounted on a movable arm 
located just above the cathode inside the liquid argon 
active volume. 

In addition, for the first time in a dual-phase detector, 
we used 37 Ar as a diffused source throughout the active 
region. 

37 Ar 

was obtained from neutron irradiation of 
nat Ar at McClellan Nuclear Research Center, similarly 
to what described in [9]. nat Ar is composed by 40 Ar 
(99.6%), 36 Ar (0.34%) and 38 Ar (0.06%). The irradia- 
tion therefore produces mainly 41 Ar that quickly decays 
with a half-life of 110 minutes, leaving 37 Ar as the only 
significant radioactive product. 37 Ar decays by electron 
capture (EC) to the ground state of 37 C1 with a half life 
of 35 days. The atomic shell cascade processes result in a 
total energy release of 2.82 keV for capture on the K-shell 
and 0.27 keV for capture on the L-shell [10]. For our pro- 
duction, 1 liter of nat Ar at 11 bar was placed next to the 
reactor core and irradiated for 4 hours to produce 30 /iCi 
of 3 Ar. The activity was verified by measuring gamma 
emission from 41 Ar shortly after irradiation. The gas was 
then cryogenically extracted, transferred to a lecturer's 
cylinder and pressurized with nat Ar to 90 bar. The acti- 
vated argon is injected in the detector trough a purifier 
(SAES MC1500-903) in a similar way as the standard ar- 
gon. In order to reach the desired activity in the detector, 
a fixed volume of gas at known pressure was introduced 
several time. It is worth mentioning that a source based 
on 37 Ar also emits monoenergetic neutrinos of 811 keV, 
as considered for example in [TTJHT^] . In the following, 
we will focus on demonstrating the usefulness of an 37 Ar 
source for low energy calibration of electromagnetic re- 
coils in dual-phase detectors. 



III. RESULTS 

We report the results of our first measurements with 
37 Ar. These measurements were taken with -Edrift = 2.4 
kV/cm applied electric field across the active region, and 
6.0 kV/cm across the 7 mm of liquid argon in the extrac- 
tion region, which ensures ~100% transmission of elec- 
trons from the liquid in to the gas phase [13] . The electric 
field in the gaseous amplification region was 9.0 kV/cm. 

The waveform of each PMT was independently digi- 
tized and stored for off-line analysis. The trigger was 
set to acquire all four PMTs signals when two selected 
PMTs cross within 5 /is a given hardware threshold set 
just above the baseline noise. Each waveform is 100 /is 
long and includes a 35 /is pre-trigger. 

The analysis of PMTs signals is performed following a 
gate-time algorithm similar to the one outlined in [14] . 
From each triggered waveform, the baseline is first com- 
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FIG. 2. Comparison of the spectra collected before (blue line) 
and after (black line) the injection of 37 Ar in the detector. 
The rightmost peak present in both spectra comes from the 
55 Fe located inside the active volume. Data obtained with 
an electric field of 2.4 kV/cm in the drift region and of 9.0 
kV/cm in the amplification region. 



puted for each channel; software thresholds are then 
set based on the single photo-electron (s.p.e.) response 
of each PMT. All digitized points above threshold are 
grouped into a single event if their time gap is less than 
a given gate time across all 4 channels. A gate time of 
3.5 [is is used for the analysis of the data reported here. 
Varying the gate time between 2 [is and 7.5 [is affects the 
results by < 1%. The single photo-electron amplitude of 
each PMT is used to compute the total event energy in 
unit of s.p.e. across all channels. 

The PMT's collection efficiency for secondary light 
varies by 30-40% for edge events with respect to cen- 
tral ones, simply because of geometry. In order to reduce 
this and other edge effects, a simple fiducialization algo- 
rithm was used to quantify the X-Y position of any event. 
For each of the four PMTs the number of detected p.e. 
is computed. The first fiducialization parameter is ex- 
tracted from the ratio of the two PMTs on the east side 
vs the two PMTs on the west side; the second parame- 
ter from the north/south ratio. The ratios are taken so 
that the parameter is always greater than 1 and are then 
mapped to a square phase space of side 1. In the follow- 
ing a fiducialization cut that selects the central ^40% of 
the active region volume is used. 

Two additional basic cuts are used for event selection. 
First of all we require that the event start time is within 
5 [is from the trigger position. We also request that no 
more than 3 s.p.e. are present in the pre-trigger and no 
more than 5 s.p.e. after the end of the event. 

In Figj2]we show the spectra collected before and after 
injecting ~2 kBq of 37 Ar in our detector. Before injec- 
tion, only the peak from the 55 Fe source is present. After 
the 37 Ar liquefied and mixed with the rest of the argon, 
two more peaks appeared from K- and L-shell electron 
capture of 37 Ar. 

The spectrum with 37 Ar is fitted in Fig(3] The low- 
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FIG. 3. Fit to the spectrum of the data collected after in- 
jection of 37 Ar. The peaks are due to the 55 Fe and 37 Ar 
calibration sources. The total fit is shown in black and the 
individual components of the fit as red dashed lines. Data 
obtained with an electric field of 2.4 kV/cm in the drift re- 
gion and of 9.0 kV/cm in the amplification region. xVd.o.f. 
= 125.5/115. 



energy background is modeled with two exponentials. 
The 55 Fe source is described by two gaussians, one at 
5.90 keV and the second at 6.49 keV. These values come 
from the weighted average of the four most prominent 
x-ray lines in the 55 Fe decay [THj. The relative amplitude 
of the two peaks was constrained accordingly to a factor 
(2.99/25.4). The width of the smaller peak at 6.49 keV 
was also fixed from the value of the primary 5.90 keV 
peak. The two 37 Ar peaks were also modeled with two 
independent gaussians. 

A relative branching ratio of 0.116 ± 0.013 for the L- 
over K-shell EC for 37 Ar is calculated from the fit. This 
is in agreement with previous measurements in gaseous 
argon P^HTTJ], 

During the same run we were able to measure the de- 
tector response to single ionization electrons (i.e.). Single 
electron events appeared following electric discharge in 
the detector, their rate decreasing over the course of sev- 
eral hours. We are currently working to study and under- 
stand the production mechanism of this single electrons. 
Here we show in in Figj4] the single electron spectrum. 
In order to study the population of single electrons, a 
stricter cut is applied, requiring less than 1 s.p.e. before 
or after each event. The fiducialization cut was removed 
as it is not effective at this very low number of p.e.. The 
average event width for single electrons was ^6 [is, with 
longer events due to pile-up, as shown in the inset of 
Figj4] The spectrum is fitted using two gaussians to de- 
scribe single- and double-electrons events. The mean and 
width of the double-electrons gaussian are constrained 
from the values of the single electrons distribution. The 
single electron response is 8.2 p.e/i.e. with a 1 sigma 
resolution of 3.4 p.e.. 

Using the single electron response, we compute the 
number of electrons extracted from liquid argon for each 
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where a and m are the recombination and mobility co- 
efficients, E the electric field and a the box dimen- 
sion parameter. Following the suggestion in [4], the 
number of initial electrons iVj is calculated assuming a 
value u> 9 = 19.5 eV for the energy required to create a 
quanta (either ionization or excitation) in liquid argon 
[5] and an initial partitioning between the two channels 
of N e JNi = 0.21 [H] : 



FIG. 4. Single electron spectrum obtained in the same elec- 
tric field configuration of the 37 Ar data. The spectrum is 
fitted with two gaussians describing the single- and double- 
electrons events. Data obtained with an electric field of 2.4 
kV/cm in the drift region and of 9.0 kV/cm in the amplifica- 
tion region. \ 2 /d.o.f. = 115.4/94. In the inset, scatter plot 
of the event width vs the collected light in p.e. for the events 
in the spectrum. 
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FIG. 5. Dependence of detected ionization in the Ar and 
55 Fe peaks from the spectrum of Fig[3]on the induced electron 
recoil energy. The solid curve is the result of the fit of the 
experimental data using the Thomas-Imel box model of eq. 
|l]) and pi). The value of the single free fit parameter is 
provided in the figure. Data obtained at an electric field of 
2.4 kV/cm. 



of the three main peaks in Fig. [3] These are plotted in 
Fig. [5] as a function of the nominal energy deposited in 
the detector. 
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(1 + N^/Ni) 
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Equation ([!]) is fitted to the data with a/(Aa 2 uE) as a 
constant free parameter, as shown in Fig. [5] Excellent 
agreement of the model with experimental data suggests 
that it is applicable in liquid argon at low energies. An 
experimental verification of this behavior at different val- 
ues of the drift electric field is currently being performed. 

It is somewhat surprising that the data agrees with the 
model up to 6 keV. Previous work in liquid xenon |22j 
has shown that above a few keV electron recoil energy, 
recombination fluctuations begin to significantly modify 
the number of collected electrons. This is manifested 
as a departure from the Thomas-Imel prediction. Our 
data offers a preliminary suggestion that in liquid argon, 
recombination fluctuations are less significant in the few 
keV energy range. 

In conclusion, we have shown that dual-phase argon 
ionization detectors are sensitive to sub-keV electron re- 
coils. We have also provided a novel calibration technique 
in this energy range through the use of 37 Ar. This first 
demonstration of sensitivity in the sub-keV to few keV 
range in liquid argon enables the possibility of a sensitive 
search for axion-like particles, via the axio-electric effect 
[23] . We are also developing a technique to produce 37 Ar 
with no 39 Ar contamination. If this is successful, this cal- 
ibration tool should be useful for liquid xenon detectors 
such as LUX [24]. Moreover, these data indicate that 
low-energy electron recoils in liquid argon can be mod- 
eled using a simple approach to electron recombination 
based on the Thomas-Imel box model. Efforts are under- 
way to further corroborate these conclusions at different 
electric fields. We are also working to get experimental 
data on the ionization yield of low-energy nuclear recoils 
for the future development of both the CNNS and dark 
matter detectors. 



IV. DISCUSSION 
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